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l -af ADAPTIVE CONTROI. LAW RFFfNLlVIKNT 
AND SOF'rW AHi: DKV Kr.OPMFNT 


Cfury I,. liiirtiTUinn 
Guntor Stein 
IIONKYWKT.L, INC. 

SIOC'TION 1 
INTHODUCTION 


I ho NASA Dryden Flight Research Center Is currently flight testing a 
ujgital fly-liy 'Wire (DFHW) flight control system installed in an l'-8C 
airi't’afl. This serves as a test vehich* for demonstrations of advanced 
control laws ami redundancy concepts to improve the performance and/or 
overall effectiveness of future flight control. In support of atlvttnced con- 
trol law efforts, Honeywell conducted a design program to define a digital 
adaptive control law suitable for flight test. The initial study (Reference 1) 
recommended an adaptive concept which combines gain-scheduled control 
laws with explicit maximum likelihood identification to provide the sched- 
uling variables. This approach was selected from a comparison of three 
camiidatc concepts: 

• Implicit gain adjustment based on self-excited limit cycles, 

• Gain adjustment based on explicit identification using a 
Liapunov model tracker, and 


• CUiin ad.juHlmeni on itlontifitaUon with Maximum 

I .ikolitioovi Hk lima (Ion. 

I .ator tIcHign t'xtj>nsionH (Hoforom t* 2) atkied a two* level esilnsate of gust 
inlenHilv and pr’ovided a new parameter update method based on Kalman 
fslimalion t»f lime- varying parameters. 

This study provided further development of the Ftirallel ('hannel Maximum 
I ik» lihot)tl Kstimation (I’CMI.K) tleslgn. A number of features have V»een 
adiled to faeilitate flight U'sting of the algorithm. The software was 
(irigmallv designed for on-board implementation. For eonvenienee and 
fit xibillty in testing, he algorithm has been implemented on the NASA/ 
DFHC Remotely Augmented Vehicle (RAV) facility (Reference 3). As 
shown in I'lgure 1, the PCMI.F software resides in a ground computer. 

The measurements requir ed by the PCMI-E algorithm- -pitch rate, normal 
acceleration, and horiz ontal stabilator position- -are received by the 

ground-based computer via the lelemeti’y downlink. The PCMLE estimates 
i!h‘ aircraft characUcn/ation parameters and compiitt'S a dynamic prt'ssurt' 
csumaic IS a liiu'ar function ot . 'I'his tpianliis is transmitted to the 
inplcN on-t>oai d digital computer and used for gain scheduling. As part of 
the groundiHilcH, mcasuia'ments wci-e restrii ted to rate gyros, accadci'o- 
meters, and servo position. Air data wei-c excliuicd l>ecaus(> aircraft like 
the F-tir, whose pc'rformancc rt'quircmenls c-an be met with air-dat;i- 
scheduled control laws, benefit most froiti adaptive control through the 
elimination of air-dutti schedules. 

The next section contains the list of symbols used throughout th'S report. 
Section 3 describes the PCMLF algorithm. The implementtition of the 
algorithm and its acceptance test are summarized in Section 4. In 
















Section 5, sensor noise statistics are developed from ground test and flight 
test data. Section 6 contains the parameter estimates from the PCMLE 
algorithm computed from recorded F-8C flight data. Cross-checks are 
provided by parameter estimates from a batch maximum likelihood algo- 
rithm. Section 7 presents recommendations '‘or flight evaluation of the 
PCMUi: algorithm. Conclusions from this study are given in Section 8. 
Three appendices to this report contain time histories of the flight maneu- 
vers and PCMLE outputs. 



SEC I TON 2 


Operators 



A 

V 


V 


2 




Superseripts 

n 

(■) 

(>“> 

( )' 


Subscripts 

( ) 


m 


( ). 


SYMPOLS 

Laplace operator 
Time derivative 
Increment 

Gradient vector with re^|A?ct to para- 
meter vector C 

Second partial derivative matrix with 
respect to parameters C 

p-th component of V(. ) 

Summation 

Absolute value 


Estimated value 
One-stop predicted value 
Value for parallel channel i 
Nominal value 


Measured value 
Value at time t 

k 


Uppcir Case* Symbols 
A 
B 

\ C»i‘ 

D 

11 

I 

.1 

1< 

\. 

L 

w 

IVl 

1V1 . M , M IM ), IVI 

q cv 5 ^ a 

M 



Diacrete system dynamics matrix 
Discrete system input matrix 
Response variable 
Measurement matrix, y due t^/ u 
Measurement matrix 
Identity matrix 

Partial likelihood function L - l/2't^lndetB 

Kalman filter gains 

Ijikelihood function 

V(>rtical gust field scale length 

Number of parallel channels 

Pitching moment coe'fficients 
indicated variables 

Trim pitching moment 

M^ value? for rigid airframe (without 
quasi -static flexibility) 

True? value of M 

f)0 

Number of data samples 
Normal and lateral acceleration 
Kalman filter covariance matrix 
A priori parameter covariance matrix 
Residual covariance matrix 


V 

V 


I'O 



Z ,Z 


vj 


0 


Sequence of N control inputs 

True air speed (m/s) 

Crossover velocity in €"• response 

Sea'icncc of N measurements 

N n ' force coefficients due to indi- 
cate' variables 


Z V 
a 


Normal force coefficient in normal 
acceleration equation 


Lower Case Symbols 

C, C. 

> 

c 

-o 

^t 

(I 

g 

h 

i 

n 

P 

q 

r 

t 

u 


Parameter vector with components c, 

A priori estimate of c 
True value of c 

Sensor displacement from c.g. (4.62m) 

Gravity 

Altitude 

Index of the minimum-L channel 
System order 
Roll rate 

Dynamic pressure 

(1) Yaw rale 

(2) Number of measurements 
Time 

Control input vector 


X 


Slate vector 


Greek Symt 
Upper Case 


Lower Case 


ft (ft), ft,.. 


Measurement vector 

Discrete system noise input matrix 

Angle-of-attack 

Gust angle-of-atlack 

Angle of sideslip 

Aerodynamic surface positions 

Dummy argument for values of parameter 
vector 0 

White noise process 

Kalman filter residuals 

White noise process 

Standard deviation of variable x 

Time constant 

Roll attitude 

Natural frequency 


SIOC riON 3 


TH’MI.K SOFTW AHi: l)K VEIX^PMKNT 


Phis Hoction summarizes the development of the Parallel Channel Maximum 
Likelihood Estimation (PCMLE) software. 3'his software implements an 
adapMve gain sehedule for the F-8C aircraft based on explicit parameter 
estimation. The software is designed for flight research using the 
Remotely Augmented Vehicle (RAV) facility at NASA/DFRC. It represents a 
further development atid refinement of an adaptive design recommended 
in References I and 2. 'Phe theoretical background relevant to this study 
IS eoniaiiu'il in these I'eports. 

A complete documentation of the PCMLE softw;ire including program 
listings anil flowcharts is available as a separate volume ^Reference 4). 

THE HASiC ALCiCtRlTHM 

The PCMUO algorithm is based on standard Maximum Likelihood Estimation 
theory as applied to longitudinal short-period F-BC dynamics. Instead of 
using the usual iterative calculations to maximize likelihood functions, 
however, it uses the parallel channel implementation shown in Figure 2. 
Several Kalman filter channels operate at fixed locations in parameter 
space. Likelihood functions are computed for each. Sensitivity equations 
are then solved only for the maximum likelihood channel and used to 
interpolate from there to the final parameter estimate with a single 
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Figure 2. Basic PCMLE Algorithm 
















Newton/ Haphnon parameter eorreetion. 'I’hls fixed strueture avoidn 
real-tinu' iterations and eliminates eor.vorgence problems. 

Theoretical identtfiability results were used to determine the number of 
parameters that could be identified with small test inputs. This accuracy 
analysis also provides insight into the number and location of the filter 
channels, 

NomiiuUly, five parallel channels are used to handle the F-8C aircraft over 
its entire operational flight envelope. The locations of these channels in 

parameter space are shown in Figure !l. Up to four paramelers-- 
surfacc effectiveness (M ), pitching moment due to ancle-of-attack (M ), 

f, 

airspeed (V), and normal force due to angle-of-attack (/ VI- -can be 

n 

estimated. I'stimation accuracy depends strongly on the signal levels in 
the control loop. I'or the small test signals producing less than 0.05 g 
HMS of normal acceU'rat ion, errors are 10 to 20 percent in M and 

-''O 

20 to 30 percent in and V which arc typical in six-degree-of-freedom 
simulation runs, 'rheoretical accuracy analyses confirm these error 
levels. 

The gain adjustment in 'he pitch and lateral control la\^^s is done on the 
basis of estimated only using scheduling functions defined in 
Hefcrence 1. However, the Ml.K design war selected in large part for its 
potential to identify additional parameters which may be needed for 
scheduling in other applications. 





Figure 3. F-8C Identifier Channel Locations 



PC’MI-K SOF'rWARE FEATURES 


The I’CMIJ' software implements the above identification algorithm in a 
flexible manner suitable for flight test experimentation. Specific options 
and features of the implementation include the following: 

• Variable number and locations of Kalman filter channels. 

• Variable Kalman filter update rates (sample skipping options). 

• Var iable number of identified parameters (up to four). 

• Output variables to monitor identification validity. 

• Provisions for four uplink parameters with fail-safe integrity 
tests, 

• Gust and rigid-body angle-of-attack estimation, 

• An optional Kalman parameter correction (instead of Newton- 
Raphson steps) for improved tracking, 

• Optional automatic adjustment of channel gains with gust 
intensity, and 

• An optional second Newton -Raphson parameter correction step. 

These features and options are achieved through a combination of software 
structure, channel model structure, and algorithm modifications. The 
software structure is discussed in detail in Section 4. Model structure 
and algorithm modifications are described below. 
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Nominal Charmul Models 


Four-slate filler models are used for each channel with the following 
discrete-time form; 

A A * 

x^=Ax +Bu + K(y - H x ) (1) 

n+l n n “^n n 

Flements of matrices A, B, K, and H are computed in a non-real time 
initialization mode from the continuous F-8C model for a specified sample 
time. 

The continuous model is 

X = I'x + Gu + r K 

( 2 ) 

>k ^ "\ * 

Individual terms appear as 




• a 

q 


I 

o 

o 

o 

i 


ai m 

q 

m 

c 

N 

z 

> a 

k 

dM dM -Z V 0 dM -Z,V 
q a a 6 6 

a a 


Op 

"b 

'e 

m m 

+ ( 


0 





(3) 


where 

M . M . M, 
q m J 

(/. V). (Z V ) 

a f 

V 
K 
d 
L 


w 


0 . 


* pitching moment derivatives 
' normal force derivatives 
true airspeed 

actuator bandwidth (12. 5 rad/sec) 
distance of accelerometer aft of c. g. 
gust field scale length (580 m) 
standard deviations of noise processes 


functions of para- 
meters to be 
identified 


Subroutine MODEL computes the discrete matrices A and H froui standard 
formulas; 


A(t) 


X“^|sl - l*i'^ 


evaluated at 
t = At 


At 

H - j A(t - t)G d T 
o 

'I'he resulting matrices have the following form: 


A 


V V V V 

V V V V 
0 0 C 0 
0 0 0 C 


p 


V 

V 
0 
C 
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where V donotew elements which vary with flight condition and C denotes 
conHtant elementH. Strictly speaking, A<3, 3)*e varies with flight 

condition also because i* depends on the bandwidth of a first-order gust 
model. However, since the incoming measurements are high-passed 
(Figure 2), the high-pass frequency, ‘“np, dominates and is therefore used 
to replace V/L . 

'I’he Kalman filler gains K and the residual covariance matrix R are defined 
by the discrete Riccatt equation; 

R (HPir^ ^ wwS 

K = (AlMl‘+ r.^wSn"’ 

rj, b'i) 

P = (A - KH IP (A - Ml ) 4 (Fjj - KW)(rj^» KW)^ 

These equations are solved via a doubly iterative algorithm carried out in 

subi'outine's DIAK and C\l . The discrete noise input matrix, Tjy is 
computed from and 1*' according to the second-order approximation 

- M I (1 + M F/2) (6) 

1) c 

In order to compute state and residual sensitivities in real time, each 

nomintil model also includes sensitivities (first partials with respect to 

parameters At , C„, C.,, C.) of the matrices A, B, H, and K. These 
o O i O 4 

are computed by numerical finite differencing techniques. 

All model parameters are initialized in non-real time and are stored in 
labelled arrays for later real-time use. 
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Mgorilhm MoiUficationH 


iVI.HUficiiiion of the banic PCMLE algorithm of Figure 2 were required to 
implement the following three options provided by the PCMLE software; 

• Kalman parameter eorrections, 

• Automatic gain adjustment with gus< intensity, and 

• rwo-step Newton- Ruphson parameters corrections. 


The first two of these options were studied under the Design Extension 
Study (Reference 2) of the original PCMLE design program. They led to 
the algorithm modifications summarized in Figures 4 and 5. Further 
details can be found in Reference 2. The third option implements a 
iv.'o-siep parameter cori’cclion by introducing a "roving" channel 
located at point c which is the estimate' obtained from the first Newton- 
Raphson step. Likelihood functions and sensitivities are computed for 
the rttving citanncl and art' used to provide the second upciate step: 


^ A 2 , , 

c = c - ( V r.l V L 


e B c 


These algorithm modifications are immarized in Figure 6. 


Models for the Second Newton- Raphson Step 

Because the second Newton-Raphson step involves a channel which is not fixed 
in parameter space, the filters and sensitivity models for this channel must be 
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Figure 4. Identifier with Kalman Filter Parameter Corrections 



















KCUMULATICK 


I 



Figure 6. Identifier with Two-Step Parameter Correction 






updated in real time. To make this process computationally efficient, the 
following approximations of the fixed channel modeling procedure are used 
in the software; 

• Matrices A and B are approximated as described below. 

• Gain matrices K a»e approximated from current min-L 
channel data; 

K = K(i ) +Vk| . , (c - c^ ) 

• Sensitivities for the approximated A, B and for H are explicit. 

• Sensitivities for K are taken from the current min-L channel; 

VK =Vk|,, 

. A i '1= 

These approximations are reasonable for small A = c - c and should not 
adversely affect the ability of the second Newton-Raphson step to estimate 
parameters. If A becomes very large, of course, the approximations 
deteriorate. It is even possible to gel u(\stable A and K combinations. 

To protect against this possibility, the second Newton-Raphson step is 
automatically bypassed if its likelihood function diverges. 

The A and B matrices used for the second Newton-Raphson step have the 
following approximated form; 



MqA^ 


V 

(^A 

A = 

V 

A~tZ^ 

V 

^A 


0 

0 

c 

0 


0 

0 

0 
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B = 
whore 

A = -y <1 + exp (-12. 5 At) ) 

B = (1 - exp (-12.5 At; ) 

M and Z are corrected for flexibility 

The seven terms circled in the A and B matrices are updated in real time. 
Elements A (1,1), (1,2), and A (2,2) use a first-order approximation in At. 
Elements A (1,4), A (2,4), B (1, 1), and B (2, 1) have essentially a second- 
order approximation. The variation of three terms in the A matrix, 

A (1, A (2, 1), A (2,3), is neglected. These are small terms with 
small percentage variations and are fixed at the model values for which 
the second Newton-Raphson step was initialized. 

In order to provide flexibility for tuning the above model approximations, 
the coefficients M^, M^, IVl Z^, Z^V can be adjusted independently from 

the nominal channel values. Their values are set through the common 
block F8 MODE which is read via NAMELIST VARL at initialization. This 
feature was specifically used to modify the function at supersonic 
conditions to better match damping of the exact and approximated models. 
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SECTION 4 


PCMLE SOFTWARE IMPI.EMENTATION 


Phis section tloscril)es the software structure, software management and 
control procedures, and acceptance tost procedures which were used to 
implement the Pt'MLE algorithm on DFRC's RAV facility. 

The PCMLE software was developed for and successfully verified on the 
Control Data CYBER 73-28 computer at NASA/DFRC. All the software 
is written in standard FORTRAN IV ;md is intended to be transferrable 
to the RAV computer (Varian V-73) for eventual flight test experiments. 

SOF'PWARE STRUCTURE 

An overview of the PCMLE software organization is shown in Figure 7. 
The computations arc divided into a background (non-real-time) segment 
to define and initialize Kalman filter channels and a real-time segment 
to process sensor data for parameter identification. Calculations per- 
formed in each of these segments are divided among a number of sub- 
I’outines, as listed in Table 1. The functions of each segment and their 
input/output structures are briefly described below. The core required 
for PCMLE is 5655 locations for subroutines plus 2730 locations for 
storage arrays. 
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INITIALIZATION 


C ENTER A 
NRTIC J 

~~r 


REAL-TIME 


READ INPUT 
DATA 


COMPUTE 

CHANNEL MODELS 
AND 

SENSITIVITIES 


INITIALIZE 
FILTERS, 
GRADIENTS. 
SKQND EARTiAL& 


WRITE OUT 
CHANNEL DATA 






C ENTER A 
PCMLE I 

IHIGH-PASS MEASUREMENTS | 

r ^ 

kALMAN FILTER FOR EACH CHANNELI 

* 

OMPUTE SENSITIVITIES FOR MIN-CHANNEL 


(OPTION) 


SENSITIVITY 
ACCUMULATION 
FOR NEWTON-RAPHSON 
PARAMETER CORRECTIONS 


1 


MEASUREMENT 
ACCUMULATION 
FOR KALMAN FILTER 
PARAMETER CORRECTIONS! 


(LOW RATE COMPUTATIONS) 


SUBCYCLE 


1. MIN-L SELECTION 

2. SIGNIFICANCE TEST AND CHANNEL CHANGE LOGIC 

3. CHANNEL DATA TRANSFER 

4. PARAMETER INCREMENT CALCULATION (NEWTON- 
RAPHSON STEP OR KALMAN FILTER UPDATE 

5. PARAMETER CORRECTION 

6. MODEL UPDATE FOR SECOND loPTlONl 

NEWTON-RAPHSON STEP vuriiur.; 

7. SECOND NEWTON-RAPHSON STEP 


4 " 

^ RETURN ^ 


Figure 7. PCMLE Software Structure 



TABLE' 1. PCMLE SUBROUTINES 


Non-Real-Time 

Subroutines 

Functions 

Core 

Required 

(decimal) 

NR PK' 

Main executive routine for non-real- 
time operation. Reads data to define 
number and location of channels, num- 
ber of parameters estimated, sample 
rate, etc. Performs all initialization 
with calls to other subroutines. 

930 

MODEL 

Defines the system matrices and sen- 
sitivities for the discrete four-state 
model described in Section 2. 

638 

FHir 

Computes high-pass filter coefficients 
and initializes filter states for each 
measurement to be high-passed. 

13 

DIAK 

CAI. 

Solves Ricatti equations for the Kalman 
filter gains of a discrete system, using 
double iteration procedures. 

101 

POLES. QRCAIT. 
QR, HESSEN 

Computes eigenvalues for channel 
models and their Kalman filter 
dynamics. 

730 

Real-Time 

Subroutines 

PCMI>E 

Main executive routine for parallel chan- 
nel MLE real-time computations. 

1767 

FH 

High-pass filter applied to measurements. 

14 

TSIG 

Produces test signal and two random 
numbers for simulated sensor noise. 

41 

FILT 

Performs fourth order Kalman filter 
update computation. 

95 
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TABLE 1. - Concluded. 


Real -'Pi me 
Subroutines 

Functions 

Core 

Requirec 

(decimal 

SENS 

Performs a sensitivity filter update for 
a given parameter. 

176 

ACCNR 

Accumulates likelihood gradients and 
approximate second partials for a 
Newton-Raphson parameter correction. 

177 

SENS2 

Performs sensitivity filter updates for 
"roving" channel of second Newton- 
Raphson step. 

176 

ACCNR2 

Accumulates likelihood gradients and 
approximate second partials for second 
Newton-Raphson step. 

177 

ACCK 

Accumulates measurements for a Kalman 
filter parameter correction. 

32 

KBF 

Performs a Kalman filter parameter 
correction. 

288 


Initialization 


The initialization of PCMLE is performed in non-real time with a call to 
subroutine NRTIC. This subroutine reads the input data deck and user 
options (UX and LX arrays) and checks the input data for reasonableness. 
It then defines the specified numbers of channels, each at its specified 
parameter values. Each channel is a four-state Kalman filter. The 
states are pitch rate, total angle-of-attack, gust angle-of-attack, and 
elevator surface position. The two measurements are pitch rate and 



normal acceleration. The input to the filter is elevator servo position. 

Two sets of gains are computed and stored for each channel, corresponding 
to low and high turbulence levels. Sensitivities are computed for each 
channel and gust level by individually perturbing each of four parameters 
to be estimated. Eigenvalues are computed for each channel model and 
each Kalman filter. All computations are performed for the sample rate 
specified. 

Real-Time Operation 

All real-time compulations are executed with CALL PCMLE (IT, lOUT), 
where lOUT is a four-component vector for uplink. lOUT(l) is the uplink 
parameter used for gain scheduling. The scalar IT is not currently used 
but would be available to further partition the real-time calculations 
should this be required by the RAV executive. PCMLE is called once per 
sample time. On the CYBER 73-28 computer, each call required 
3 to 6 msec of computer time. 

During each call, the sampled values of pitch rate, normal acceleration, and 
elevator servo position are high-passed. Residuals and likelihood functions 
are computed for each channel (fixed in parameter space). Gradients 
and second partials are accumulated if Newton-Raphson parameter 
corrections are selected. Otherwise, measurements are defined for a 
Kalman filter parameter correction. The remaining real-time operations 
are spread over seven subcycles executed in sequence, as shown in 
P’igure 7. During real-time operation, the detailed performance of the 
algoriit m can be monitored by the UX array. This is defined later in the 
subsection on outputs from PCMLE,, 
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If the sample skipping option is selected, the sample time is modified 
uttomatically by NRTIC, However, PCMLE must then be called at the 
appropriate new sample rate. Either one or two samples may be skipped 
in this manner. 

Real-Time Inputs to PCMLE > -Real-time inputs to PCMLE from the RAV 
program are assigned to user array elements UX(ll) through UX(18) and 
UX(20). l)X(ll) through UX(13) hold the usual measurements needed by 
the identifier; UX(14) and UX(15) are used to compute an average elevator 
servo position when the sample skipping option is used. UX(16), UX(17) 

and UX(18) are used to communicate "true" values of M , , M , and 

6o a 

q when these quantities are supplied by a simulation, and, finally, UX(20) 
provides real-time adjustment of test signal magnitude. These inputs 
and their units are summarized in Table 2. 

Real-Time Outputs from PCMLE - -Provision has been made for four 

outputs to be supplied via the calling argument lOlIT. At present only 

JOlJ'l’(l) is used. It is a 1 /q estimated from M. and scaled to be 

he 

between 0 and 512. (Scale factor is 50000. ) 

Other outputs from PCMLE are contained in the real-time IJX array. The 
scaled test signal is in UX(19). The remaining outputs are in UX(1) through 
IIX(IO) and IJX(21) through UX(50). These are primarily used for 
monitoring PCMLE performance. All PCMI.E outputs are defined in 
Table 3. 
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TAHLE 2. HEAL-TIME INPUT VARIABLE ASSIGNMENT 


User 


Mnemonic 


Variables 

Description 

Expression 

Units 

I'Ndl) 

Pitch rate 

Y (1) 

rad sec * 

12 

Normal acceleration 

Y (2) 

-2 

it sec 


Servo position at cur- 
rent lime t, 
k 

DELTA (1) 

rad 

14 

Servo position at t^^ ^ 

DEI.TA (2) 

rad 

15 

Servo position at tj^ ^ 

DELTA CD 

rad 

16 

RAV q estimate 

QUART 

lbs ft ^ 

17 

RAVM, estimate 
6o 

MDOr 

-2 

sec 

18 

RAV M estimate 

Cl 

MAT 

-2 

sec 

20 

Test signal magnitude 

SIGUT = SIGH TO 
UX (20) 

0 SIGUT 40 I f 

-2 

tt sec 












Mnemonic Expression 1 Units 


* h 

N CC 


^2 

fu 
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Several halt conditions have been defined for PCMLE. The*ie are listed in 
Table 4. The first two stops occur in the initialization mode if the data 
(loos not define a stable filter. The last three stops occur in real time if 
the parameter update option is not defined. In all the test cases run, these 
slops have never been encountered. However, in the final flight test 
softwar e, they may be switched to a mode change operation rather than 
a halt. 


TABLE 4. PCMLE HALT CONDITIONS 


Slop 

Condition 

41 

Ricatti etiuation not converging in subroutine CAL during 
initialization (Unstable Model). 

31 

Inverse does not exist in subroutine DIAK for computer filter 
gains. 

Check data deck. 

21 

No inverse exists in Newton-Raphson parameter update. 

22 

'.(o inverse exists in second Newton-Raphson parameter update. 

11 

No invers • Ts in Kalman parameter update. 


User Inputs an d (Outputs 

Both the background (initialization) and real-time program segments pro 
vide user inputs and options. These inputs are selected with a nominal 
data deck and/or by setting the UX and LX user arrays in the RAV 
software. 
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i hc' I \ . I l i as st rvt's ilifl'ei tMU tumMioriH dep<*iuiing on whelhcM' (he program 
iH III tiu' re il tiiin* or iii the initiali/alioti segment. Hoth funetions are 
defineit (x'losv. 

hipuls - - l'’ollo\ving loading, the WAV program (eontaining all (he 
IH !\ll,l'; subrouiines) exeemes the noti-real-time initiali/.ation segment. 
During the first siieh initialization operation, subroutine NHTK’ reads a 
latvi deek ami stores it on a ttanfiorarv tbsk tile for quiek restarts. 

I’hi> input jiarameters vshieh are read-in on data ixmls are did’ined in 
'I'able f>. rhe ilata dei k defiiuxs nominal (liefaiilt) vxihu's for all program 
paraim'liM's inehuimg the five ehantu'l loeations tabulateil in Tabh* (1 (and 
illustrated in l''iguri> ID. 


t'onsoh^ Inputs Prior to Iniliali/ation -Certain Uigieal variables u '^s) 
di'fiiu'd in Table 7 mas b( assigned from iht' eont rol eonsole prior to 
initialization. If this step is bvpasstni, the P('M1,I'' default option will 
be thi' basie IH’lMId'l algorithm (T'igure ID with a single Nmvton* Haphson 
parameter eorrei lion. These logieals I'aiinot bi' altert'd in real time. 'To 
ehange them ihi' usm‘ must leave real tinu' atui rmiiitializ.t'. 


It is also possible to rmlefiiu' ei'rtain data deek parameters (irior to 
inil lali/at u>n by setting U\ t'b'menls from the eonsob\ 'These eUnhents 
are listiui in Tabh' 8. 'The redefinition is performed as followst 


\ ariable 


Nominal 
default 
value from 
liata deek 


( UX arras'^ 
value J 
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TABLE 5. NAMEIJST PARAMETERS 


Fiirnmeter 


NX 


n r PRINT 

NC 

NPO 

-ISO 

IGO 

WUTO 


Definition 

Nominal 

Value 

Number of states in Kalman filter (q, oi, mg, fi) 

4. 

Number of measurements (q, Nz) 

2. 

Number of noise sources (ag, q sensor, 
Nz sensor, 5 sensor) 

4. 

Maximum number of measurements 

3. 

Convergence criteria in DIAK for Ricatti 
solution 

10 -« 

Maximum number of iterations in DIAK for 
Ricatti solution 

20. 

Nominal sample time (automatically increased 
if samples arc skipped) 

0. 02 

Print interval when LX(5) is true (available in 
batch mode only) 

0. 5 

Number of channels 

5. ! 

Number of parameters estimated in first 
Newton- Raphson step 

2. 

Starting channel location 

3. 

Starting gust level 

1 . 

Natural frequency of shaping filter on random 
test signal 

6. 

Damping of shaping filter on random test signal 

1.25 









Parameter 

Definition 

Nominal 

Value 

SIGUTO 

Gain of shaping filter on random test signal 

4. 

IITIO 

initial state of shaping filter on random tost 
signal 

0. 

UT20 

Initial state of shaping filter on random test 
signal 

(above state initialization is useful for gene- 
rating deterministic square waves or sine waves 
with TSIG subroutine) 

0. 

UTMAX 

Magnitude limit of filtered test signal 

10. 

TAUPO 

Time constant of likelihood accumulation filter 

5. 

TAUP2 

2 

Time constant of low-pass filter on L, VL, 7 L 

0. 6 

WHP 

Cutoff frequency of second-order high-pass on 
measurements 

2. 

DHP 

Damping of second-order high -pass on 
measurements 

0. 7 

NP20 

Number of parameters in second NR update 

2. 

SIGGO 

RMS statistic assumed for gyro noise (deg/scc) 

0. 15 

SIGACCO 

RMS statistic assumed for accelerometer 
noise (g) 

0.02 

SIGWLO 

RMS statistic assumed for low w gust level 
(ft/sec) 

1.0 

SIGWHO 

RMS statistic assumed for high w gust level 
(ft/sec) 

5.0 







TABl.E 5. - Continued 




Nominal 

Parameter 

Definition 

Value 

nvsTo 

HMS statistic assumed for elevator servo 
hysteresis (rad) 

lO""* 

n 

'I’hrcshold parameters controlling gust level 

0.01 

'mrrrjGO ( 

switch 

3.22 

irr.KX') 

Threshold parameters controlling channel 

0. 1 

riitrrjct^j 

switch 

3. 22 

rrrjzo 

threshold parameters controlling ZIMIN 

0. 1 

rilRT.JZOJ 

sc'lcction 

13. 8 

IMMSO 

Threshold parameter for significance of 
likelihood function 

0.25 

ZPl MAX 

Maximum value limit of parameter 1 estimate 

-1. 

ZP2 MAX 

Maximum value limit of parameter 2 estimate 

1.3 

ZP3 MAX 

Maximum value limit of parameter 3 estimate 

120. 

ZP4 MAX 

Maximum value limit of parameter 4 estimate 

10. 

ZPl MIN 

Minimum value limit of parameter 1 estimate 

-75. 

ZP2 MIN 

Minimum value limit of parameter 2 estimate 

-0. 3 

ZP3 MIN 

Minimum value limit of parameter 3 estimate 

-60. 

ZP4 MIN 

Minimum value limit of parameter 4 estimate 

-10. 
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TABLE 5. - Continued 


Parameter 

Definition 

Nominal 

Value 

GSCiLO 

2 

Diagonal elements of V 1^ matrix (four numbers) 

0.001 
0. 1 

0.0003 
0. 01 

ZP 

Matrix defining location of channels (4 para- 

See 

Table fi 


meters per channel x 10 channels max = 
40 numbers) 

NHKF 

Number of measurements in KF parameter 
I'orrection (q, N/.) 

2. 

NXKl' 

Number of states in KF parameter model 
(Zj.Z.,,Z.^,Z4,Zi) 

5. 

NIIKF 

Nuntber of gradient measurements available 
in KF parameter correction 

4. 

rlKl’O ^ 

Parameters in KF slate model 

0. 001 

T2KI-' J 


1.0 

NCYC 

Number of subcycles 

i 

7. 

giK") 

Initial covariance matrix for KF parameter 

-6 


correction (five numbers) 

5 X 10 
0. 2 

0.0045 
0. 00a 

ISKPO 

Number of measurements skipped (0, 1 or 2) 

0. 


NSO 


Signul-to-noise ratio (increased residual noise 
statistics at high signal levels) 


0. 05 






PAHLE 5. - ConcJudod 


Parameter 


Definition 


Nominal 

Value 


IVIQO 


!\1Q12 


M \ ! '2 


Coefficients which define channel models from 
C.,, C.^. via; 

- MQO + (MQl + IVIQ12 C2) 

M = (M AI 4 MAI 2 : (’2) M, 

ft 6 

V (VI + C2) ]/ 

b 

Z V (ZAVl 4 C4) M. 


-0.23 


0.03 


0. 61 


0. 92 


Z \V1 


Z ZDVl 'M 


ZDVl 


DIST 


Distance Nz is measured aft of c.g. (ft) 


15. 15 


AC'mW 


ZPIC' 


Paiuiwidth of first-order actuator model for 
^ (rad) 

Quasi-static flexibility corrections ; 


M, - (1. 4 M (FXl 4 M, FX2)) 

^ b o ho 

ifM <ZP1C M, - (1. 4 ZPIC (FXl 

ho h 


4 (ZPIC) FX2 


12. 5 


Scale length in gust model (ft) 


1750. 


QME 


Scale factor between q (psf) and M 

(q - QME M, ) ^ 

ho 


ALPHWL 


Bias between water line a and zero lift n (rad) 


•M 



TABLE 6. NOMINAI. TIIANNEl, LOCA I’lON 


( h.intu'l 


Pafami'ti't’s 

o 


^’3 

‘4 

-2. 34 

0 

0 

1 

0 

-5. 27 

0 

0 

0 

- n . !> 

0 

0 

0 

-26. 7 

0 

0 

0 

-26,7 

i 

i. 

60. 

0 


l ABLE 7. TNITIALIZA TION USEB LOC.ICAL ASSK'.NMENTS 


llsi't- 
1 .oRU'al 

I'lou'tion 

Mnomonu' 

LX(1) 

Not UH«'(I 


2 

Noi ustai 


3 

KnabU's Kalm:m pafanu'tc'i' rorriTiion in 
plant' t)l' lii'sl NovMon- llaphson I'ufri'ft ion 

KAl. 

4 

Not us I'll 


5 

\ol US.'.I 


6 

Knablt's stn'oiui Nnwton- Raph.son slop 

NR 2 

7 

Not UHt'tl 


a 

Inhibits initial nuiili'l printout during 
mitiali/ation 

Noiu' 

10 

Atitls psouilo- ranilt)m sonsor noi.sf to 
J'CMl K'.s input incasuri'im'ntK 

NOS S 

a -50 

Not ust'tl 



Di'laull ValiU's 











TABLE 8. INITIALIZATION USER VARIABLE ASSIGNMENTS 



Sign . 







FARI.E 8. - Concluded 









ShouUi the atIjusU'd variable Tall outside of program imposed limits, the 
variable is returned to the default value. It is also important to note that 
RA\ resets all I’X array values to zero whenever the program leaves 
real time, itent e, desired changes of the nominal data deck must be 
reentered into the l)X array before each reinitialization. 

C onsole Inputs du ring Real- Time Operation - -The logicals defined in 
Table 9 can be set and reset during real-time operation. If they are set 
in non- real time, then the corresponding mneumonie will be defined in 
the PCM 1, 10 subroutine as soon as real time is entered. 

Only one user variable, IIX<20), can bo set from the console in real time. 
Its function is to alter PCMI.lO’s test signal magnitude (see Table 2). 

SOPTWARK CONFIGURA TION CONTROL PLAN 

lOffective on the date of acceptance at DFRC, a software configuration 
comrol plan was implemented to protect the integrity of the PCMLE 
program. The control plan includes the following elements: 

• Source file management 

• Change reporting and execution 

• Verification test 


Each of these elements is discussed below. 


PABLE 9. REAL-TIME USER LCX3ICAL ASSIGNMENTS 


User 

Logicals 

Function 

Mnemonic 

LX(l) 

Disables automatic channel changes 

NCX'IIC 

2 

Disables first Newton-Raphson or Kalman 
stop 

NOEST 

‘A 

Not used 


4 

Disables automatic gust level changes 

NOCIIG 

5 

Enables printout from PCMLE (use in 
batch mode only) 

PRINT 

6 

Not used 


7 

Not used 


8 

Not used 


9 

Enables PCMLE algorithm 

MLE 

10-50 

Not used 



Default Values; LX(I) 








iCMi’ Manage ment 

FUt?s of controlled software documentation were established following veri- 
fication of program operation on the RAV facility. Duplicate files are 
being maintained at both DFRC and Honeywell, Each file contains; 

• PCMLE source deck 

• Flowcharts 

• Listings (current plus two previous versions) 

• User information 

• Change notices 

To maintain a record of the software, changes in any item require a 
properly executed change notice. 

ha nge Reporting 

The change notice required for controlled file alteration is a standard form 
having the following parts; 

A) Reason for change- -Brief summary of problem, 

B) Description of change- -Specific program changes, documentation 
changes, and affected procedures. Mark-up copies of documents 
affected by change are attached. 

C) Functional checkout- -Provisions of checkout procedures re- 
quired to verify correct operation of modification. 
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D) Program change report --Listing of who performed the changes, 
data perforn.od, and verification procedures followed. 

Each change notice must be signed off by the originator and approved by 
Honeywell’s program manager, A sample change request form is shown 
in Figure 8. 

V'erification Test 

Verification tests are those procedures that are performed to assure that 
software changes are accomplished in the intended manner. The following 
actions are taken: 

• Following change in the master source deck, a now listing is run 
and checked against the change notice requirements, 

• For critical changes, portions of the acceptance test routines 
are run. These requirements are listed in part (C) of the 
Change Request, 

PCMLE ACCEPTANCE TEST 


This section summarizes the procedures followed in conducting the PCMLE 
Acceptance Test on the F-8C simulation at NASA Dryden Flight Research 
Center, Test conditions, selected time histories, and interpretation are 
given below. 
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Program: PCMLE 


FORTRAN SOFTWARE PROGRAM CHANGE REQUEST 


Originator 


Date 


Approval 


Effectivity 


Title of Change 


A) Reason for Change 


B) Description ot Change 


Change 
Statement No. 


Statement 


C) Functional Checkout 


D) Change Report 


Who performed 


Change /Checkout 


Date 


Figure 8, Software Change Request Form 













Dofinition of Tost Cases 

Aceeptaiiee test cast’s were jointly defined by DI*'HC and Honeywell 
pefsonnel, I'hey are summarized in Vuble 10, 

I )e tailed I'e s t P roeedur e and C o nditio ns 

Tile acceptance test was run on the CDC Cyber 73-28 Computer in real time. 
The l' -8C simulation with the control laws function us one program, the 
HAV executive with PC'MlvH as another. Data interchange was accomplished 
via 1)/A and A/l) trunk liiu’s as presently mt’Clumized in the RAV mode. 


rije acceptance tests used a standard acceleration maneuver. The l’'-8C 
simulation was brought to the (5100 m. 250 KlAS flight condition and trimmed 
to 1 g flight, riu' aircraft was ttien accelerated at constant altitude until 
it stabilized at a new velocity (approximately 300 m/sec). The simulation 
was run man-in-tlie loop using the l'’-8C Iron Hird. Small pilot inputs 
wa't'i' us('d as required to maintain trim. 


Two i’ighi -clumnei strip-cliart recorders were used to obtain time histories 
of variinis aircraft response variables and selectml PC!MLK outputs, 
riiese figures are presented later in this section, Por the performance 
evaluation we have tin? luxury of knowing "true” values of and M , 
rtu'se parameters were computed in real time by approximati'lv deter- 
mining tlu' slope of the and Cjvi simulation functions. These approxi- 
malt'd slopt's w’ere ustni with PCMLE estimates to compute the and 
1\1 errors shown in the time histories. 
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TAB*.E 10. TEST CASES FOR ACCEPTANCE TEST 


'l\‘St 

Case 

Description 

1 

Baseline --Identification of M C2. agust» a rigid body 

of the simulated F-8C. The maneuver was a full after- 
burner acceleration from 250 KIAS to 450 KIAS at 6100 m 
in the presence of 2 m/sec gust disturbance and no sensor 
noisi% Five channels located as defined in Table 6 were 
used. The nominal sample rate was 50 sps. 

2 

Repetition of baseline in absence of gusts to show auto- 
matic adjustment of Kalman filter gains. 

3 

Repetition of Case 2 with a smaller number of parallel 
channels (4), 

4 

Repetition of Case 2 with a slower Kalman filter update rate 
(two samples skipped). 

5 

Repetition of Case 2 with Kalman filter parameter correc- 
tions instead of Newton -Raphson corrections. 

(i 

Rep(‘tition of Case 2 with two Newton- Raphson iteration steps. 

7 

Repetition of Case 2 with two additional parameters 
estimated (C^ and C^), 


The detailc?d test procedure is given in Table 11 and the acceptance 
criteria are shown in Table 12. 


A cc ep tance Test Res ults 
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A successful acceptance test was completed at DFRC on 15 October 1978. 
Table 13 summarizes the time histories for each of the test cases. In 




some instances several additional time histories are included to illustrate 
the performance of the algorithm. 

The performance of the algorithm is judged by examining the q error 
traces and the and error traces. For most cases the performance 
is as expected and compares with results in Reference 2. For the sample 
skipping cases there is an increase in the errors when the high gust 
channel is used. This effect warrants further investigation. 


TABLE 11. ACCEPTANCE TEST PROCEDURE 


Test 

Case 

Procedure 

1 

PCMLE was initialized by reading the nominal input data cards 
(Reference 4). The program was then put into real-time oper- 
ation. After 10 seconds, estimation was armed via logical 
MLE [l.X(9), Table 9J and the accelerated maneuver was 
initiated. Variables appropriate to judging PCMLE were 
monitored in real time via CRT displays of UX and LX arrays 
and via strip eliart recordings. 

2 

Same procedure as Baseline Test Case 1 is followed with no 
turbulence in the F-8C simulation. 

3 

PCMLE was reinitialized prior to this run. Only four channels 
w<mt^ used, Channel 1 was deleted and Channels 2 through 5 
were renumb(?red as 1 tlu'ough 4. This verifies that fewer num- 
bers of channels can be used and that their location can be 
moved. Thesii changes were accomplished with a new input 
data deck. Following initialization, PCMLE was brought into 
real time and engaged as discussed under the Baseline case. 

4 

The I^CMLE program was reinitialized with the sample skipping 
optional 1PX(28) = 2, Table 8j. Every third measurement of q 
and Nz m(?asur(?ment was used, and the 6^. servo position 
measurements W(>re used to compute an average value valid over 
the thrt!e-samplt? interval (Reference 4). The sample rate of the 
Kalman filters was one third of the simulation rate. Real-time 
operation was the same as discussed under the Baseline case?. 

5 

PCMLE was reinitialized with the logical KAL set (LX(3), 
Table 7j and then brought into real time as discussed under 
the Baseline case. 
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TABLE 11, - Concluded 


Procedure 


PCMLE was reinitialized with the NR2 logical set lLX(6), 

Table 7j and the NOEST logical lLX(2), Table 9J. This 
combination holds the first Newton -Raphson stop at the chan- 
nel location, and, hence, the second Newton-Raphson step 
should approximate the first stop of the Baseline case, PCMLE 
was then brought into real time as discussed above. Two para- 
meters were estimated by the second step. 


PCMLE was reinitialized with NP = 4 lUX(32) = 2, Table 8j 
to estimate two additional parameters. Real-time operation 
was the same as the baseline. 



TABLE 12. ACCEPTANCE CRITERIA 


Criteria 

Qualitative --The channel switching shall be smooth and in the 
proper sequence (channels ordered in more negative 
The parameter estimates shall be smooth during channel 
switches. 

Quantitative --The error in estimated dynamic pressure shall 
b«‘ within ± 50 percent of the " true" value throughout the 
transition. 

The channel indicator shall show a switch to low gains when 
data p(?rmits, although it is not required to remain continuously 
on the low gust channel. Qualitative characteristics of base- 
line apply. 

Sinci' th(' first channel was not selected as a min-L channel 
in the baseline run, the rc'sults will match the baseline run 
except for renumbering of the channels (i. e. , 1 was 2 in base- 
line, 2 was 5, etc.) 

Only qualitativi' factors of baselim? apply. Errors will approxi- 
mately match baseline. 

This run shall demonstrate improved tracking over the base- 
lint' case (less error in estimated dynamic pressure). 

'I’he estimates will approximately follow the baseline run. 
(Differences are due to model approximations included in the 
second Newton-Raphson computations.) 

This run shall produce estimates of C and with the 
qualitative characteristics of the baseline. 






















Figure 8A, Aircraft Response --Test Case 1 
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Figure 8B. PCMLE Response --Test Case 1 
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Figure 9A. Aircraft Response --Test Case 2 
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Figure 9B, PCMLE Response --Test Case 2 
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Figure IID. PCMLE Response --Test Case 4 with Sensor Noise 















Figure 12B. PCMLE ilespon.so--Test Case 5 
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Figure 14B. PC'MLE Hospomse 


Test Case 6 




Figure 14C. Aircraft R(‘sponse- -Test Case fi Deceleration 
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Figure 14D. PCMLE Response- -Test Case (> Deceleration 


SKCriON 5 


SKNSOH NOISl*: MODKUNC. 


I'hf PC’MLK algorithm was di'Vi'lopeU under the assumption that sensor 
tioisc stati.sLies are eonstanf over the flight envelope ami arc reasonably 
well known. Hence, they are not trcjited as parameters to be identified by 
the algorithm. I’hc valiiiity of tins assuTnption is investigated in this 
seeiion. Sensor tiata from grouttd tests at engine off, idle, and 80 percent 
maximum HPM, and from flight tests are analyzed. Ht'suUs support the 
assumption that the statistics are constant, but modified nominal values 
are ixHjuired to match the tf'st aircr.il’t's effective sensor i’haractt?ristics. 

KNCilN'K-Ol’P DATA 

An 82-secomi segment of sampled gvro and acceleroim'lt'r outputs under 
ciuiesccnl hatigar conditions was analv/eii. I'ht* analysc's includi'd mean 
ami variance calculations, histogram plots, and powcu’ specnral dcnsitii's 
(PSDs) cotnputed via Past Fouriei' Transform methods. The variance 
calcadaiions arc' summai’t/ed m I'able 14. I’he> show RMS h'vels 
roughl> ('c(ua) to one-third of each sen,->(M''s least -significant cjuanti- 
/aijon hi! (I.Slll. I'hc corrcsponiling histograms arc plotted in Figures la 
and 18 and the I’SDs arc given in Figures 17 and 18. These show that 
noise in the hangar is dominated by relatively white random motions of 
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Direct transformation of 4090 data points, with resulting plots smoothed 
by averaging adjacent frequency samples. 


I’ABUO 14. S I ANDARD DEVIATION OF SENSOR NOISE 



Ground Test Data I 

Sensor 

Engine Off 

Engine Idle 

80% RPM 

.Normal accelerometer 
1 2 -bit resolution 
I. .SB ^ 0. 00547 g 

0.0018 

0.0042 

0.016 

J’ilch rate gyro 
12 -bit resolution 
I. SB = 0. 038 deg /see 

0.011 

0.019 

0.057 


the last one or two quantization bits. Fsin^ a Gaussian tissumption for the 
underlying noise processes which move these bits, RMS levels of roughly 
one-third bit are again obtained for both sensors. Since there is this much 
similarity between the gyro and accelerometer noise levels, it appears that 
the bit motions arc generated by A/D electronics rather than by internally 
generated sensor noise. 

ENGINE-ON DATA 

I’he RMS sensor outputs increase substantially when the engine is running. 
'I'his is shown in Table 14 for two engine spceds--idle and 80 percent RPM. 
PSDs for these conditions arc shown in Figures 19 through 22. They indi- 
cate that most of the RMS increase can be traced directly to various 
resonances between 3 and 20 Hz. While these resonances may in fact be 
legitimate input signals as far as the instruments are concerned (i.e. , not 
internal sensor noise), they must be treated as "effective sensor noise" 
for purposes of PCMLE because the algorithm includes no models to explain 
the sensed motion. 
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Figure 15. Gyro Noise Distribution (Engine Off) 
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Figure 16. Accelerometer Noise Distribution (Engine Off) 
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Figure 20. Accelerometer PSD (Engine I 
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FLIGHT DATA 


Sorisor noise numbers applicable once the aircraft leaves the ground were 
ded iced from Kalman filter residual histories, as generated during flight 
data parameter identification runs discussed in Section 6. These para- 
meter identification runs used a general purpose identification program 
(GPMLE) to fit a best linear model to the flight data. The model extracts 
estimates of the (rigid body) sensor output, leaving the residuals, 

A 

” ■'^k* either internal sensor noise, unmodeled dynamics 

(e.g. , structural modes), atmospheric turbulence, and mismatched rigid 
body motion. In the aoscnce of turbulence and assuming a good model fit 
for rigid body motion, therefore, the residuals provide "effective" sensor 
noise time histories directly. 

HMS noise levels from such residual histories are summarized in Table 15. 
Four maneuvers are shown, corresponding to flight data segments docu- 
mented in Section 6. It is evident from this table that effective noise 
numbers in flight are substantially higher than both the hangar data and 
the ground test data. This is highlighted in Figures 23 and 24 which illus- 
trate the data from all test conditions in graphical form. The figures 
clearly show that ground and hangar tests are inadequate indicators of 
airborne noise statistics. They also show that, while the constant statistics 
assumption made for PCMLE seems reasonably valid, the actual RMS 
levels used in the acceptance test should be modified somewhat to match 
the test aircraft sensors. 


83 


r \m,K 15. S rAM)\R|) dkviation ok kkkkc tiv k skxsor 

\01SK; Kl.ICIlir DATA 


Vlaneuvt't* 

ti :2 

2:1 

2:4 

2:6 

Altitude (fet't) 

20200 

20129 

19H4B 

22500 

Mach 

0.44 

0. 566 

0. 85 

1.12 

GPMl.K Residual 
1 evels (RMSi 





Aci-elerometer 

(g’st 

0. 05tt 

0. 050 

0.050 

0.054 

C'lV ro (deg /sec) 

0. 11 

0.079 

0. 12 

0.13 


As a iui'thi'f evaluation of effective noise statistics, it would have been use 
ful to I’epeat the analyst's in 'I’able 15 untler turbulence conditions. There 
is some rationale to suggest lh;it effective sensor noise should increase 
further with turbulence level because of increased unmodeled structural 
t'xcitation. In the abst'nce of turbulence data, we are forced to rely on 
conclusions from relatetl noise modeling efforts conducted on the SAAH 
,I A-57 aircraft (Reference 5). These suggest that effective noise incretises 
due to turbulence* are probable negligible. 














SECTION 6 


PARAMETER ESTIMATION WITH FLIGHT DATA 


In order to gain inereused confidence in the algorithm prior to flight test, 
PCMLE WiTs exercised with recorded sensor outputs from F-8C flight 
tests. Since the flight recordings do not contain PCMLE's own test signal, 
data segments with large pilot commands were used to provide good 
conditions for identification. Results of these off-line exercises are very 
positive and provide a high level of confidence for successVul closed-loop 
flight tests. 

In addition to the PCMLE exercises, a general purpose maximum likelihood 
estimation (GPMLE) algorithm (ieveloped for the F-8C (Reference 1) was 
used to estimate all pitch axis paramoters in a conventional iterative batch- 
processing mode. Results from both the PCMLE and GPMLE estimation 
are presented and compared in this section. They provide a data base for 
flight test recommendations made later in the report. 

TEST POINTS 


Flight conditions for which flight data were processed are plotted in Figure 
25. This is an adequate number of conditions for checking PCMLE, although 
more data at 40,000 ft (12,195 M) would have been desirable. The flight 
records examined contain 15 pitch doublet maneuvers covering a dynamic 
range from 126 psf to 840 psf. These have all been processed with PCMLE. 
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Figure 25. F-8C Flight Data Test Points 



Several lateral maneuvers were also used as "disturbances" to PCMLE to 
evaluate the effect of lateral maneuvers on the pitch axis estimation. 
Finally, an acceleration run from Mach = 0. 82 at 37,000 ft to Mach = 1. 15 
at 30, 000 ft was processed. 

The flight data for each test point consist of time histories for the three 
measurements needed to drive the PCMLE software. The time histories 
were sampled at 50 sps. Other related measurement parameters are 
given in Table 16. Note that the accelerometer was located at the c.g. 
Also, the quantization level of each sensor is higher than the ground data 
quantization used in Section 5. 
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TABLE 16. F-8C MEASUREMENTS 


Measurement 

Prefilter 

Scale 

Quantization 

Pilch rate (q) 

40 llz 

+ 70 deg/sec 

0. 138 deg/sec 

Normal acceleration 

40 Hz 

±8 g 

0.0156 g 

Elevator position 

'V 

40 Hz 

+11.75% 

-21.5^ 

0. 0325 deg 


Six longitudinal maneuvers (two consecutive pilot doublets each) were 
extracted from the Flight H2 data tape. Each maneuver is a 10. 24-second 
segment. The maneuver start times are given in Table 17 for the time 
references of the tape. Five similar segments containing pitch axis pilot 
commands were extracted from the Flight #3 data tape. These are also 
identified in Table 17. For shorthand reference to all i2 data segments, 
the symbol "Maneuver i;j" will be used to designate Flight i, Maneuver j. 

Plots of the maneuver time histories for all 11 maneuvers are shown in 
F’igures A1 through All in Appendix A. Note that the accelerometer 
measurement on Maneuvers 2:5 and 3;5 are contaminated with low frequency 
oscillations. 

PCMLE PERFORMANCE 

The above maneuvers were used to exercise the baseline PCMLE algorithm 
(two parameters were identified using a single Newton-Raphson step) and 
also two of the software 


options : 
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TABLE 17. FLIGHT TEST POINTS 




1. The second Newton-Raphson correction, and 

2. Identification of additional parameters. 

Since the algorithm for all these eases is nominally set up to operate with 
measured servo position, a modified actuator bandwidth parameter, ACTBW, 
was used to operate with the available measured surface position instead. The 
bandwidth was moved from 12.5 rad/sec to 100 rad/sec. In addition, the accel- 
erometer noise pai'ameter, SIGACCO. was increased from the nominal value of 
0.02 g RMS to 0.04 g RMS, as suggested in Section 5. 'Ibe distance param- 
eter. DIST, was set to zero to match the c.g. location of the instrument. 

All other parameters remained at their nominal values in Table 5. 

The Kalman filters in PCMI.K use a parameterization based on "unflexed" M^ 
because the functions are simpler. Therefore the parameters estimated by 
fitting data to the model will also be "unflexed. " 

Table 18 summarizes performance of the baseline algorithm for the 11 
flight test points. Estimated variables in the table are the following: 

A 

“ PCMLE's estimate of surface effectiveness before quasi - 
state flexibility corrections. 

A 

- Surface effectiveness estimate after quasi -state flexibility 
corrections. 

- PCMLE’s estimate of the one-sigma accuracy of its M 5 ^ 

^ estimate. This tends to be an optimistic number because 

PCMLE does not recognize errors due to unidentified 
parameters. 

A 

- Estimate of small perturbation parameter C., used to calcu- 
late pitching moment due to angle-of-attack. 
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TABLE 18. PERFORMANCE OF PCMLE ON FLIGHT DATA 




















V 

CTl 

Gust 

I.ovel 

SIGSQ 


Estimate of the one- sigma accuracy of the estimate. 

Scaled up by 1M5^I, this parameter gives a rough (optimistic) 
indication of the expected accuracy of PCMLE's estimate. 

Estimate of pitching moment coefficient due to angle-of- 

A ^ 

attack, as computed from M5 and C„, 

Estimate of velocity 

PCMLE's selected min-L channel. 

PCMLE's estimate of random vertical gust level, either higi? 
(5 ft/sec RMS) or low (1 ft /sec RMS). 

Estimated scale factor on the residual magnitudes of the 
min-L channel. SIGSQ = 1 corresponds to nominal noise 
conditions. 


Time histories of the min-L channel's gyro and accelerometer residuals, 

A A A A A 

the estimated RMS error •"'ud the Mj^, M5^^, M^^ and V estimates 

themselves are shown for each maneuver in Figures B1 through Bll of 
Appendix B. Note that the starting transients include some drift in the 
estimates since PCMLE is not getting any information until pilot commands 
start. As mentioned earlier, the normal PCMLE test signal is not 
present in any of the flight (iata. 

Compared with expected (simulation) parameter values for the test points 
in Table 17, all the estimates in Table 18 are reasonable except those 

A 

for Maneuver 2;5. This case produces a more negative M^^ estimate than 
expected, especially when compared to Maneuver 2:6 which is nearly the 
same flight condition. Looking at the raw data for Maneuver 2:5, we see 
that the accelerometer is particularly noisy for this maneuver and contains 
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the unexplained low Iroqueney OHCillntions noted earlier (about 3,25 Hz). 
This probably explains why PCMLE goes to the high gust estimate and 

selects a lower Mr value. Note, however, that the amount of shift in Mr 

"o '^o 

(50 percent over the value on Maneuver 2:6) should cause no closed-loop 
stability or performance problems. 

Following the above baseline r*in3, a selected subset of cases was rerun 
with a single channel located at the parameter estimates from the first run. 
rhis corresponds to a second Newton-Raphson correction performed in 
sequential fashion. Results of these experiments are summarized in Table 
19. Their time histories are shown in Appendix C. 

The results verify two properties of the PCMLE algorithm; 

1. The first Newton-Raphson parameter correction achieves 

imp'^oved fit to the flight data. This is evident by comparing 
residual traces for the baseline cases with residual traces from 
the second iteration. The baseline residuals correspond to the 
min-L channels indicated in Table 18, while the second residuals 
correspond to channels located at corrected parameter values 
from the first Newton-Raphson steps. Note that the second 
residuals are smaller but still do not resemble w/iite noise 
during the pilot input periods. This is because PCMLE's channel 
models ignore several aircraft parameters which are weakly 
identifiable under test signal conditions but can produce sub- 
stantial residual errors under large pilot inputs. 
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2. The first Newton- Raphson correction comes close to achieving a 
local minimum of the two-parameter likelihood function imple- 
mented in baseline PCMLE. This is evidenced by the fact that 
the second corrections in Table 19 do not move far from the first 
corrections in Table 18. Most changes are within one- or two- 
sigma units of the algorithm's own optimistic accuracy estimates. 
Hence, the second step is "satisfied" with the location found by 
the first. We note again that, because PCMl.E's models ignore 
several other aircraft parameters, the location of this local 
minimum does not necessarily correspond to the true parameter 
values. According to our identifiability and design studies, how- 
ever, it should be accurate to within 10 percent or so. 

As a final exp<‘riment, the effects pf estimating additional parameters were 
examined using Maneuvers 2;8 and 2;5. Results of these tests are summar- 
ized in Table 20. For each maneuver, two-, three-, and four-parameter 
identification trials were run. These show small changes ("dative to 
PCMI.E's accuracy estimates) of the original two parameters when addi- 
tional parameters are estimated. The additional parameters themselves 
are found only crudely, as indicated by their corresponding accuracy esti- 
mates. F’or example, the expected one-sigma error on C., (small pertur- 

iJ 

bation paramct<>r for velocity) is greater than 35. The maximum variations 
of C.^ are known from wind tunnel data to be only +60. Similarly, the small 
perturbation parameter forZ^V has expected one-sigma errors greater than 
7. 20. Its maximum variations are known to be + 10. 0. Hence, while 
PCMLE produces numbers for the additional parameters, their accuracy 
is hardly better than a priori knowledge. This is consistent with past 
identifiability and design studies. 
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TABLE 20. VARYING NUMBER OF PARAMETERS ESTIMATED (One Channel Located 
at Corrected Parameter Values from First IteratitMi) 



























PARAMETER ESITMATION WITH GPMLE 


Solcctod maneuvors from Table 18 were also used to conduct identification 
runs with a general purpose maximum likelihood algorithm (GPMLE). This 
software was used to do the original identifiability analyses which led to the 
PCMLE design (Reference 1). It uses a conventional iterative batch- 
processing approach to parameter estimation. The parameter estimates 
arc updated with standard Newton- Raphson steps until the likelihood function 
ceases to improve. For this algorithm, the various gradients required are 
analytically computed. 


Ident ific a tion Model s for GPMLE 


The identification model used by GPMLE was a three-state pitch axis model 
with discrete measurements of pitch I’ute and normal acceleration. The 
model is 
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where 


? , T] = white noise 

d = acceleration displacement from c.g. 


Modtl Parumetcmalion 


One of the feutures of the PCMLK algorithm that allows it to work well 
while estimating a small number of parameters is its method of parameter- 
i/.ution. The coefficients appearing in the above model are computed from 
one dominant parameter (which is M plus iither small pei turbation 
parameters (C through C., C’„> as shown in Tabie 21. 


TABLE 21. I'-aC Mf)UEl. PARAMETERIZATION 


-0. 22 + (0. 028 - 0. on Cg) + C j 

(0. 61 + 0. 92 C.J C. 

& a 


(200 + C^)f^ 


Z V = (53 + CO 
a 4 5 


Z V “ (7. 7 + C-) C. 
ho 6 5 


Correction for quasi-static flexibility: 

M. = M. (1 +0.016 M. + 0.0002 Mt) 

he 0 o ho ho 

ZfV = (Z.V ) M. /M, 

h 6 o he 0 o 





OPMU*: Results 


The parameter estimates obtained with GPMLE are summarized in Table 22. 
Kstimates oi‘ M^, V, and Z^V are plotted against in Figures 

26 through 20, The figures are the original scatter plots used in the model 
analysis on the F-8C adaptive study (Reference 1). The x’s and o's repre- 
sent the model parameter values at 25 flight conditions which were used to 
establish the PCMLE functions given in Table 21. The functions are plotted 
as solid lines. The four flight data points from Table 22 are plotted as A's. 

Comparison of the flight data with the original linear models shows that the 
model fits quite well for V, and Z^V (Figures 27, 28, and 30), How- 
ever, the aircraft seems to have more damping (M^, Figure 26) and larger 
c.g. acceleration due to surface deflection (ZgV, FTgure 30) than indicated 
by the model. Moreover, there does not appear to be a need for the quasi- 
static flexibility correction used in Table 21. The original model data in 
Figures 27 through 30 are plotted as a function of "unflexed” surface effec- 
tiveness, M^^, while the llight data is plotted as a function of actual 
("flexed") 

Since these plots are compatible (except for the scale factor changes on M 
and already mentioned), it follows that mav just as well be 
interpreted as M 5 . 
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Figure 28. 







GPMLE and PCMLK Compar isons 

A comparison of PCMLE and GPMLI3 estimates for a subset of maneuvers 
is given in Table 23. For the and parameters, a percent difference 
is also shown which indicates that the two estimation procedures agree to 
within reasonable percentages. True values of these parameters are, of 
course not known. For the V estimate, on the other hand, a measured air 
ilata value is available. This value falls somewhere between the GPMLE and 
PCMI.E estimates. The GPMl-E velocity estimate does not improve over the 
PCMI.E estimate at two test points <2:2 and 3:1), Both the percentage differ- 
ences and the differem'es between estimated and measured values are consis- 
tent with theoretical performance predicted during the PCMLE design program. 

Hopresentative comparisons of GPMLE and PCMLE residual time histories 
(v - y-y) are shown in Figure 31. These traces correspond to the residuals 
from the last iteri tion of GPMLE (a Kalman filter located at the parameter 
values in 'Pable 22) as compared with the residuals from the min-L channel 
of PCMi.E (a Kalman filter located at one of the nominal channel locations 
given in Table 6). Both filters fit the raw signals quite well. However, it 
is clear that the GPMLE filter should (and does) fit better because it 
includes several aircraft model parameters not recognized by the PCMLE 
filters. The net effects of this improved fit are the 10 to 20 percent 
parameter differences already noted in Table 23. 

MANEUVERING FLIGHT 

The performance of PCMLE during a maneuver is shown in Figure 32. The 
top five traces show the response of the aircraft. The maneuver, lasting 
about 135 seconds, is an acceleration from Mach = 0. 85 to Mach = 1.15 
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TABLE 23. GPMLE/PCMLE COMPARISONS 
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Figure 31. GPMLF./PCMLE Residual Comparisons (concluded) 








Figure 32. PCMLE Performance in Maneuvering Flight 



during which the altitude decreases from 41,000 feet to 23,000 feet. This 
is immediately followed by a deceleration back to Mach « 0. 89. The F- 8 C 
is supersonic for about 57 seconds during this maneuver. 

A A 

The next two trac’es show the Ma and estimates from PCMLE. Note 

A ® ® 

how goes sharply more negative (as it should) as the aircraft goes super- 

A 

“onic, ITie M 5 estimate was used to produce an estimated dynamic pressure 
q . In the bottom trace of Figure 32, this estimate is compared to a dynamic 

pressure (q ) computed from the measured altitude and mach number. 

n 

The q e rror is initially large (for 10 seconds or so) because there is no 
pilot activity. (This maneuver does not contain any test signal. ) During 
the remainder of the maneuver the RMS error is about 20 percent. 






SECTION 7 


FLIGHT TEST RECOMMENDATIONS 


Based on the acceptance test and flight data processing results discussed 
in previous sections, the PCMLE software is judged to be ready for flight 
test evaluation. Recommended nominal parameters, test inputs, and 
(‘Valuation experiments which should be incorporated in the flight tests are 
discussed in this section. 

NOMINAL PCMLE PARAMETERS 


The flight tests should be initiated with the same nominal PCMLE para- 
meters recorded in Table 5, except for the following modifications; 


1. Modified PCMLE model parameterization 


MQl - 0. 044 
ZDVl =13.8 
FXl = 0. 
FX2 = 0. 


(old value 0.028) 
(old value 7, 7) 
(old value 0. 016) 
(old value 0. u002) 


The first two changes alter the pitch damping function and the acceleration 
due to surface deflection function in PCMLE's models. The remaining 
changes remove quasi-static flexibility. These are justified by GPMLE 
results in Section 6. 
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2 . 


Modified sensor noise 


SIGACCO ■ 0.06 


(old value 0. 02) 


This change increases the accelerometer noise level to the value found in 


TEST SIGNALS 

Test inputs for PCMLE flight evaluation should include normal pilot inputs 
ap.i flight disturbances, plus random test signals generated by the filter 
network shown in Figure 33. 


RANDOf! 

NUMBER 

r.ENERATOR 



LIMITER 


Figure 33. Test Signal Generation 


This filter is mechanized as a subroutine of PCMLE and corresponds to the 
same test signal routine used for simulator design evaluation at Langley 
Research Center. The output should be applied as a C*-command to the 
pitch axis control augmentation system. The random number generator 
provides a uniform distribution from -0. 5 to +0. 5. The frequency and 
damping of the filter are adjustable, as is the RMS level of the test signal. 
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Provisions have been added to produce sine wave and square wave test 
signals with the same software. To do this the input and damping are set 
to zero and initial conditions on the filter are defined to produce a sine 
wave. By limiting the sine wave, a suitable square wave can be realized. 

A FOR '^1 AN coded version of the test signal subroutine is shown in Table 
24. The various parameters are communicated through labeled common 
(UTEST). Typical values are given in Table 25, 


TABLE 24. FORTH/ N CODE FOR TEST SIGNAL 


SUBRO'iTlNF TSin(MOr»Ef'»T) 

COMMON/tlTEST/nT*SFfn*WUT*OUT*‘‘ IGUr*iJTin#OT20tUTMAX 

C UT/INPtjT s S*<;TGuT/(S#‘> ♦ ?«0(iT«ttuT*S ♦ «ijT*wun 

C TMRUT B RAMOOM NO W[Tm 0 MFaN UNfFO*»M -0,5 TO 0,5 

C IIT1,UT? STATFS in SFrONO 

Of) TO ^•^) Mooe 
C TNITIAI IZATION 

1 W2t)Ts-WUT«W»JT*nT 
T?WUTa-?,*WUT*OUT*nT 
OAMUTOsOOWT (-?A.*T7WOT) 
tlTlsUTlO 

UT?=ur?o 

return 

2 5FFDsAmoD(31?5.»SFF0*14 35n71«1<S5. ) 
5sO,2R103H30A5fr7F-lO*SFFI> - 0.5 
5sW?ur*UTl ♦ T7wi)T*tlT? ♦ OA vUT0»S 10 JT*S 
UTlsUTl ♦ nT»(»T? 

UT?sUT? ♦ S 
UTsUT? 

IF (AR'?(ilT) ,OF,MTMAv) ilT = SIOM (MTMAX , JT > 

return 

ENO 
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TABLE 25. TEST SIGNAL PARAMETERS 


«) (WUT) » 6,0 --Frequency 
C (DUT) = 1.25 --Damping ratio 
SIGUT «4.0 --RMS level 
Initial Values: 


SEED 

= 3051758125 

--Random number generator seed 

IITIO 

= 0 

--Initial test signal 

UT20 

= 0 

--Initial test signal rate 

IJTMAX 

= 10 

--Test signal magnitude limit 


FLIGHT EXPERIMENTS 

Recommended flight test experiments with the PCMLE software are sum- 
marized in Table 26. They fall into seven major groups which should be 
completed in sequential fashion in order to maximize safety and experi- 
mental value. An eighth group of experiments which involves off-line 
processing of data from other groups is also recommended. These off- 
line runs can be conducted as data become available. They serve to maxi- 
mize the experimental value of available flight hours. A brief description 
of each experimental group is given in the following pages. 

f 

I' 

I 

f 
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ccptabili ly 


These experiments ean be condueted with the standard air-data scheduled 
C -command control mode. The random test signal should be mechanized 
as shown above and inserted at the C ' -command point. With several differ- 
ent pilots and at several straight and level trimmed conditions, the RMS 
test signal level should be slowly increased by a ground-based experimenter 
from zero to 10 ft /sec. The pilot should be asked to indicate when he first 
detects the signal and when it becomes unacceptable. The corresponding 
levels should be noted by the experimenter. We recognize that the defini- 
tions of detectability and acceptability must necessarily remain vague and 
that the entire experiment will at best be *' informal" in the human factors 
sense. More sophisticated experimentation is not justified unless the test 
signal level tiu ns out to be a crucial design and performance issue. 


G rou p 2 ; Open- Loop RA V Opera t icm 


These experiments repeat selected maneuvers from this report in the real- 
time KAV' environment but with uplink data unused (i. e. , control gains 
should be set by the air-data schedule). Each experiment should quali- 
tatively match estimation results presented in this report. This will verify 
proper real-time downlink/RAV/PCMLE operation. The received uplink 
gain parameters (unused) should be compared via telemetry with on-board 
air-data scheduled gains and with PCMLE's sent uplink parameters. This 
verifies proper uplink and gain-scheduled operation. 
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Group 3; Open- L oop RAV Operatio n With Tes t Signals 


These experiments examine estimation accuracy in straight and level "hands 
off" flight at several test signal levels. The aircraft should be in the stan- 
dard air-diUa scheduled C-Command Mode, and PCMLE should be in its 
baseline configuration. The test signal level should range from less than 
delectable to barely acceptable, as determined from Group 1 experiments. 

Group 4: Clos ed “L oop RAV Operation 

These experiments repeat selected maneuvers from this report with the 
RAV loop closed. That is. the aircraft should be in the Adaptive C ^-Command 
Mode, and PCMLE should be in Baseline. Cases with pilot commands only 
and with pilot commands plus selected test signal levels (from Group 3) 
should be run and should qualitatively match estimation results presented 
in this report. Closed-loop handling qualities should be judged by the 
pilot and should closely approximate the scheduled C*CAS mode ratings. 

If these flight results are positive, other test points not covered in this 
report should be evaluated as available. 

G rou p 5; Plight Transiti ons 

These experiments examine PCMLE*s tracking ipability in closed-loop 
RAV operation. Cases shoidd be run with test inputs only and with occas- 
ional (normal) pilot inputs. 
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Group 6; All«Auitudc Maneuver in g Flight 


These experiments evaluate PCMLE performance during various flight 
maneuvers and configuration changes. The aircraft should remain in the 
Adaptive C‘i'-Commund Mode throughout, and PCMLE should be in Baseline 
with a test signal level judged acceptable from previous flights. 

Group 7;__Flight in Turbulent Ai r 

PCMLE performance should be evaluated in turbulent flight environments 
as available. 

G r oup 8 1 Off-Line Data Proc es s ing 

rhese experiments use selec ted flight data from Groups 1 through 7 lo eval- 
uate various PCMLE options. Using prerecorded data for this purpose 
serves two functions. First, it makes more effective use of available flight 
hours, and, second, it provides a bi.seline run over the same data against 
which to make performance comparisons. 

The options which show greatest promise include; 

• Channel Reconfiguration 

More or fewer channels 

More channels without Newton-Raphson parameter corrections 

• Kalman Parameter Corrections 

Unaided (as presently mechanized) 

Aided with other data (as discussed in Reference 2) 
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• Sampli? Skipping 


Indications to date are that the second Newton -Raphson step option and the 
opUon to find additional parameters have little promise and hence should be 
of lower priority. Two recommended channel reconfiguration options are 
given in Table 27. One uses only three channels and hence relies heavily 
on the Newton-Raphson parameter correction step for accurate estimation. 
The other uses ten channels. It should be evaluated with and without a 
Newton-Raphson parameter correction step and also in combination with 
sample skipping. Both procedures would offset the increased computing 
time needed to handle the large number of channels. 

Additional details on the recommended experiments in each of the above 
groups can be found in Table 26. We note that these experiments are tiot 
intended to represent a rigid pi'otocol for flight experimentation with PCMLE; 
rather, they should be viewed as a rough experimental outline to be enhanced 
and modified as the opportunities of the moment permit. 
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SECTION 8 


CONCLUSIONS 


Thc' goal of this program was to refine the PCMLE design and prepare it 
for flight evaluation on DFRC’s Remotely Augmented Vehicle facility. 

The refinement includes several steps: 

1. The addition of a number of options to enhance the research value 
of flight tests by permitting easy modification to the baseline 
configuration <such as adding channels, varying the sample rate, 
etc. ). Proper operation of these features has been verified on 
the F-8C Iron Bird. 

2. The determination of sensor noise statistics from ground teats 
and in-flight recordings. Results show that the statistics are 
reasonably constant over the flight envelope. Specific values 
recommended for the PCMLE algorithm are based on time histor- 
ies of Kalman filter residuals from the flight records , The 
assumption that sensor noise does not have to be identified on- 
line was confirmed. 

3. Identification performance was checked using the flight data. 

These estimates were cross-checked with batch MLE identi- 
fication and compared with wind tunnel data. The estimates are 
consistent. Overall performance on flight data correlates well 
with theoretical predictions and simulation results. 
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'rho PC’MLE design is now ready for flight test. Specific experiments are 
recommended in Section 7. Successful flight demonstration of the design 
signifies renewed vitality of adaptive flight controls in modern digital 
implementations . 
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APPENDIX A 


FLIGHT DATA TIME HISTORIES 
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APPENDIX B 


PCMLE PERFORMANCE TIME HISTORIES; 
BASELINE AJ XJORITIIM 
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Figure B5, PCMLE Performance, Maneuver 2:5 (concluded) 
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Conversion Factors to SI Units 


To correct from- 

To- 

Multiply by- 

ft 

rn 

0.3048 

ft /sec 

m/sec 

0.3048 

ft/sec^ 

tn/sec 

0.3048 

psf 

N/m^ 

47.88 
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